Determining scalar field potential in power-law cosmology with observational data 
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In power-law cosmology, we determine potential function of a canonical scalar field in FLRW 
universe in presence of barotropic perfect fluid. The combined WMAP5+BAO-I-SN dataset and 
WMAP5 dataset are used here to determine the value of the potential. The datasets suggest slightly 
closed universe. If the universe is closed, the exponents of the power-law cosmology are q = 1.01 
(WMAP5 dataset) and q = 0.985 (combined dataset). The lower limits of ao (closed geometry) are 
5.1 X 10^® for WMAP5 dataset and 9.85 x 10^'^ for the combined dataset. The domination of the 
power-law term over the curvature and barotropic density terms is characterised by the inflection 
of the potential curve. This happens when the universe is 5.3 Gyr old for both datasets. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The presence of a scalar field is motivated by many 
ideas in high energy physics and quantum gravities, al- 
though it has not been discovered experimentally. TeV- 
scale experiments at LHC and Tevatron may be able to 
confirm its existence. It is nevertheless widely accepted 
in several theoretical modeling frameworks, especially in 
contemporary cosmology, in which an early-time acceler- 
ated expansion, i.e., inflation, is proposed to be driven 
by a scalar field in order to solve horizon and flatness 
problems ^J. After inflation, components of barotropic 
fluids such as radiation and other non-relativistic matter 
were produced during reheating and cooling-down pro- 
cesses. A scalar fleld was also believed to be responsible 
for the present acceleration in various models of dark en- 
ergy P]. The present acceleration is strongly backed up 
by various observations, e.g. the cosmic microwave back- 
ground [3] , large-scale structure surveys |4] and SNe type 
la observations [H El El • 

Power-law cosmology, where a o^f , describes an accel- 
eration phase li q > 1. Modelling the present expansion 
with a power-law function where g ~ 1 was found to be 
consistent with nucleosynthesis O [9], the age of high- 
redshift objects such as globular clusters [5J [TUJ HH [H], 
SNe la data [U, SNe la with H{z) data [13], and X- 
ray gas mass fraction measurement of galaxy clusters 
[m [15]. In the context of the power-law model, other 
aspects such as gravitational lensing statistics [T5] , angu- 
lar size-redshift data of compact radio sources [H], and 
SN la magnitude-redshift relation [TT] [T7] have also been 
studied. Originally, the power-law expansion has its mo- 
tivation from the simplest inflationary model that can 
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remove the flatness and horizon problems with simple 
spectrum [18 . For the present universe, the idea of lin- 
ear coasting cosmology (a oc t) [I^l can resolve the age 
problem of the CDM model [TU] while as well agreeing 
with the nucleosynthesis constraint. The coasting model 
arises from non-minimally coupled scalar-tensor theory 
in which the scalar field couples to the curvature to con- 
tribute to the energy density that cancels out the vacuum 
energy [9l [20]. The model could also be a result of the 
domination of an SU(2) cosmological instanton [5T]. 

Here our assumption is that the universe is expanding 
in the form of the power law function. Two major in- 
gredients are scalar field dark energy evolving under the 
scalar field potential V{(f), and barotropic fiuid consist- 
ing of cold dark matter and baryons. We derive the po- 
tential, and use the combined WMAP5 data f22] as well 
as the WMAP5 data alone to determine the values of q 
and other relevant parameters of the potential. The nu- 
merical results are subsequently compared and discussed. 

II. COSMOLOGICAL SYSTEM WITH 
POWER-LAW EXPANSION 

Two perfect fluids, the cold dark matter and scalar 
field (f) = (t){t), in the late FLRW universe of the sim- 
plest CDM model with zero cosmological constant are 
considered. The time evolution of the barotropic fiuid 
is governed by the fiuid equation p-^ — —3Hp^, with a 
solution 



where n = 3(1 + w-y) and I? > is a proportional con- 
stant. For the scalar field, supposed that it is mini- 
mally coupled to gravity, its Lagrangian density is £ = 
(jy^ /2 — V{4>). The energy density and pressure are 

P0= + p^=\4>^ -v{<j>). (2) 
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The fluid equation of the field describing its energy con- 
servation as the universe expands is 



3H<j) 



V = 0. 



(3) 



Total energy density ptot and total pressure ptot of the 
mixture are simply the sums of those contributed by each 
fluid, for which the Fricdmann equation is 



H' = 



8ttG 



-Ptot 



tc 
7a 



It is straightforward to show that 



V{(t>) 




(4) 



(5) 



where SttG is related to the reduced Planck mass Afp by 
SttG i\/p ^. The power-law scale factor is 



a{t) = flo 



(6) 



without flxing oq = 1 at the present time because we have 
implicitly rescaled it to allow for k taking only either one 
of the three discrete values 0, ±1. The Hubble parameter 
is 



H{t) = 



a{t) 



(7) 



Our goal is to construct V{t) using recent observational 
data, as far as the simplest CDM model is concerned. 



III. SCALAR FIELD POTENTIAL 



The sign of k depends on the sign of the density pa- 
rameter f2fc = —kc^/a^H^. In our convention here, 
fc = 1 (fifc < 0) for a closed universe, fc = for a flat 
one, and k — —I (flk > 0) for an open one. The present 
value of the scale factor can be found from the definition 
of ilk,o, that is, 



ao = 



The density constant D can be found from Q, 
D = ^7,0^0 — ^f,oPc,oal, 



(11) 



(12) 



where fi^ g = ricDM.o + ^b,o, i-G- the sum of the present 
density parameters of the barotropic fiuid components. 
Pcfi is the present value of the critical density. The neu- 
trino contribution is assumed to be negligible. The values 
of Hq, to, rife^Oi f^CDM,Oi and rif,,o are taken from obser- 
vational data. 



B. Observational Data 

We work on two sets of data provided by 22J. One 
comes solely from the WMAP5 data and the other is 
the WMAP5 data combined with distance measurements 
from Type la supernovae (SN) and the Baryon Acoustic 
Oscillations (BAO) in the distribution of galaxies. For to, 
Ho, ^b,Oi and f^cDM.Oi we take their maximum likelihood 
values. The curvature density parameter ^Ik.o comes as 
a range with 95% confidence level on deviation from the 
simplest ACDM model. The data are shown in Table [l] 



We will work with observational data in SI units. 
Restoring the physical constants in place, we obtain 



H 

y 



2fcc2 
"3^ 



Dc^ 



(8) 



where Mp — Hc/SttG and we have set n — 3 {w^ — for 
dust). Incorporating (|6| and ^ into the above equation, 
we obtain 



V{t) = 



Mlc 3q^ -q 2kc% 



t^ 



alt^i 



,3q 

'•0 

a^t^'} ■ 



(9) 



We shall consider contribution of the first term alone in 
comparison to total contribution when including the sec- 
ond (the curvature) and the third (density) terms. 



A. Cosmological Parameters 

Using the equation for the Hubble parameter ([t]) at 
the present time, we have 



Hoto- 



(10) 



IV. RESULTS AND DISCUSSIONS 

Using combined WMAP5+BA0+SN dataset, the po- 
tential is 

, 1.03x1026 1.5x1023 1.5xl0'*2 
nt) = ^ + ^2s^^ (13) 

whereas, for WMAP5 dataset alone. 



V{t)^ 



1.11 X 10^ 

t2 



7.6 X 10^ 

^203 



4.6 X 10 

t^ 



43 



(14) 



in SI units. We use the mean of each ^k.o interval to rep- 
resent Clk,o in each of the above equations. Their plots are 
shown in Fig. [ij In both cases, f2fe.o is negative (a closed 
universe). The points at which the potential, its deriva- 
tive, and its second-order derivative, are zero (tintcrccpt, 
^max, and iinfloction, respectively) are also determined, for 
both ^k,o and each end of the ^k,o interval. The results 
are summarised in Table HTl 

The values of the exponent q from the two sets of data 
are only slightly different, but only the latter is an accel- 
erated expansion as g > 1. The determination of q from 
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Parameter WMAP5+BAO+SN WMAP5 

to 13.72 Gyr 13.69 Gyr 

Ho 70.2 km/s/Mpc 72.4 km/s/Mpc 

f^b.o 0.0459 0.0432 

ncDM.o 0.231 0.206 

Q.ko -0.0179 < < 0.0081 -0.063 < fi^ o < 0.017 



TABLE I: Observational data used in the construction of our scalar-field potentials |22) 





WMAP5+BAO+SN 




WMAP5 




Qufi = -0.0045 


-0.0175 < Q,kx) < 0.0085 


!^fc,o = -0.023 


-0.063 < fifc.o < 0.017 




0.985 


0.985 


1.01 


1.01 




1.9 X 10^'^ 


ao > 9.85 X 10^*^ (closed) 
ao > 1.5 X 10" (open) 


8.4 X lO^'* 


ao > 5.1 X 10^^ (closed) 
ao > 9.8 X 10^^ (open) 


^intercept 


2.7 Gyr 


2.62 Gyr < t < 2.7 Gyr 


2.7 Gyr 


2.6 Gyr < t < 2.8 Gyr 


imax 


4.0 Gyr 


3.94 Gyr < t < 4.0 Gyr 


4.0 Gyr 


3.8 Gyr < t < 4.1 Gyr 


^inflection 


5.3 Gyr 


5.26 Gyr < t < 5.4 Gyr 


5.3 Gyr 


5.1 Gyr < t < 5.5 Gyr 



TABLE II: A summary of numerical results. Times are shown in Gyr for comprehensibility. Positive and negative Qk'a 
correspond to open and closed universes, respectively. 




FIG. 1: The potentials in ([l3| and ([l4|. The units of the 
abscissa and ordinate axes are sec and J/m'^, respectively. 
The crosses mark their inflection points. Also plotted in dash 
lines are their first terms. Each potential does not actually 
converge to its first term, but later intersect with and deviate 
from it, though still very close together. However, this occurs 
much later (at t = 2.8 x lO**"* sec = 8.8 x lO''^ Gyr in both 
cases). 



X-Ray gas mass fractions in galaxy clusters favours open 
universe with q > I (q = 1.14 ± 0.05) [15 and combined 
analysis from SNLS and H{z) data (from Germini Deep 
Deep Survey) assuming open geometry yields q = 1.31 
[13j . Note that, in the power-law regime, q only depends 
on the observed values of the Hubble constant and to- 
This may give an impression that the maximum likeli- 
hood values from the combined data has yet to be relied 
upon, but the power-law expansion has not been proven 
to be the case nonetheless. 



FIG. 2: The potentials in ( |l3]l a nd (l4\ a long with the rad- 
icands of the integrands in ( |16[ ) and i\l7\ (dash line). The 
shaded region is the post-infiection phase. The unit of the 
abscissa axis is sec. After tinflcction, <?f>(i) is real. 



After ^inflection, the potential from each data behaves 
like its first term, i.e. decreasing in its value while in- 
creasing in its slope (being less and less negative). The 
other terms quickly become weaker. This can be seen 
in Fig. [l] Since the first term is contributed only by 
H(t) (and its time derivative), it is dominant in the post- 
infiection phase. In fact, the convergence to zero of the 
potential is slower than its first term alone (see ( 13 1 and 
( 14 1 ) , because the sum of the last two terms consequently 
becomes positive before converging to zero. This means 
that the plots of each potential and its first term in Fig. 
[T] eventually crosses, but it occurs much, much later at 
t = 8.8 X 10^^ Gyr. Along with the potential function in 
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^ , we also obtain the solution 



CONCLUSION 



0(t) = 



2M^c 



H - 



dt 



(15) 



in SI units. Using WMAP5+BA0+SN dataset, 



1.06 X 1026 1.5 X 1023 3.0 X 1042 



<2 tl-97 

where, for WMAP5 dataset alone, 



^2.96 



1.09 X 1026 7.6 X 1024 9.3 x 1043 



(16) 

dt. 
(17) 

In the late post-inflection phase, the first term is dom- 
inant over the k and D terms then the last two terms 
of the radicands are negligible (Fig. [2]). The above two 
equations are approximated as 



t2 



t2-l 



^3.04 



(j){t) w 1.04 X lO^lni, 
whereas, for WMAP5 dataset alone, 

(t){t) w 1.03 X 10" Int. 



(18) 



(19) 



The radicand in ([17| of the WMAP5 dataset is zero at 
approximately ^inflection — 5.3 Gyr (see Fig. [2]), therefore 
so does <t>{t). While the combined dataset has the zero 
radicand (then zero 4'{t)) in (16 1 later at approximately 
t = 5.4 Gyr. Scalar field exact solutions for the power-law 
cosmology with non-zero curvature and non-zero matter 
density are reported in [53]. It is also worth noting that 
the general exact form of the potential, that renders scal- 
ing solution, is some negative powers of a hyperbolic sine 

m. 



We consider a potential function of a homogeneous 
scalar field in late-time FLRW universe of the simplest 
CDM model with zero cosmological constant, assuming 
power-law expansion. The scalar field is minimally cou- 
pled to gravity and the other fluid is non-relativistic 
barotropic perfect fluid. We use two sets of observa- 
tional data, combined WMAP5-I-BA0-I-SN dataset and 
WMAP5 dataset, as the inputs. Potential functions are 
obtained using numerical values from the observations. 
Mean values of both sets suggest slightly closed geome- 
try. The WMAP5 dataset implies accelerated expansion 
{q = 1.01) while the combined dataset gives q = 0.985. 
This is slightly lower than the value obtained from SNLS 
and H{z) data (g = 1.31) [T3] and X-Ray gas mass frac- 
tion {q — 1.14 ± 0.05) [H]. Our result is independent of 
the geometry unlike q obtained from |13j which assumes 
open geometry. For closed universe, the WMAP5 dataset 
puts the lower limit of 5.1 x 1026 for while the com- 
bined dataset puts the lower hmit of 9.85 x 1026. We 
characterise the domination of the first term of ([9| by 
using the infiection of the potential plots from which the 
first term is found to be dominant to the potential 5.3 
Gyr after the Big Bang in both datasets. 
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